Quantum antiferromagnetism in the d = 3 Hubbard model — 
a spin-fluctuation approach 
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A self-consistent spin-fluctuation theory is developed to obtain Tn vs. U for the half-filled Hubbard 
model in the whole U /t range. Good agreement is obtained in the strong coupling limit with 
the high-temperature series-expansion result for the equivalent Heisenberg model. Quantum spin- 
fluctuation correction to the sublattice magnetization is also obtained for all U at the one-loop level. 
A spin picture is used throughout, and quantum effects are incorporated through transverse spin 
fluctuations, which are evaluated in the random phase approximation using a new method. 
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Finite-temperature antiferromagnetism within the 
three-dimensional Hubbard model at half filling involves 
two fundamentally different aspects. One is the local 
symmetry breaking at temperatures much smaller than 
the AF gap parameter A, and formation of the local- 
moment which grows rapidly and saturates with increas- 
ing U /t. The other is the scale of low-energy spin exci- 
tations, which initially rises along with the charge gap 
for small t/, and then crosses over to an approximately 

/U fall-off for large [/, becoming degenerate, in the 
strong coupling limit, with the exchange energy scale 
J = At'^/U of nearest- neighbor (NN) spin couplings in the 
equivalent S — 1/2 quantum Heisenberg model. Thus, 
while in the weak-coupling limit {U << B = 12t, the 
free-particle band-width) it is essentially the process of 
moment melting which determines the magnetic transi- 
tion temperature Tn ~ A, in the strong coupling limit 
(J7 >> -B) it is the spin disordering tendency of thermal 
fluctuations which determines the transition temperature 
Tn J at which long-range AF order is destroyed. A 
proper description of the magnetic transition in the whole 
U/t range within a single theory is therefore quite chal- 
lenging. 

Earlier studies of finite-temperature antiferromag- 
netism within the three-dimensional Hubbard model, 
aimed at determining Tn vs. U in the whole U jt range, 
have employed the functional integral formalism, mainly 
within the static approximation, quantum Monte 

Carlo methods, and recently the dynamical mean- 
field theory, Q which becomes exact in the limit of in- 
finite dimensions. All share the common feature of 
yielding a Tn which approaches, in the strong coupling 
limit, the mean-field-theory (MFT) result for the equiv- 
alent spin-1/2 Heisenberg model, T^'^ = zt'^/U, where 
z = 2d is the lattice coordination number. While 
this is natural for the functional-integral schemes em- 
ploying static approximation and the dynamical mean- 
field theory, due to neglect of the Goldstone mode, high- 
temperature series expansions yield a substantially lower 
Tn ~ iMt^/U. |,|,0§ A recent mean-field approach 
utilizing the Onsager reaction-field correction yields a Tn 



in close agreement. ||l2| Here account is made of the low- 
energy spin excitations through an approximate map- 
ping to an effective 5* = 1/2 Heisenberg model with U- 
dependent, extended-range spin couplings. 

In this paper we have considered a spin picture 
throughout for the description of quantum antiferromag- 
netism in the Hubbard model. We imagine quantum 
spins on each site, and incorporate (quantum) transverse 
spin fluctuations about the (classical) HF state, having 
local moment {Sz)'iw- Within this conceptually simple 
framework we have obtained both the Neel temperature 
Tn vs. U and the sublattice magnetization m vs. U (at 
T = 0) in the whole U/t range. Transverse spin fluc- 
tuations are obtained within the random phase approx- 
imation (RPA) in terms of magnon mode energies and 
amplitudes. A new method is described for evaluating 
the transverse spin correlations for arbitrary U jt. 

As (5*2 )hf is the maximum (classical) spin polariza- 
tion in the z-direction, and therefore also the maximum 
eigenvalue of the local Sz operator, hence S = (5z}hf 
plays the role of the effective spin quantum number. At 
T = 0, quantum corrections to the sublattice magneti- 
zation are then obtained at the one-loop level from the 
RPA- level transverse spin fluctuations. On the other 
hand, at finite temperature a self-consistent theory is 
used in which the magnon spectral functions and energies 
are self-consistently renormalized. Consistent with our 
RPA-level description in which details beyond this level 
are ignored, only momentum-independent multiplicative 
renormalizations are considered. While assumption is 
made here regarding their form, these renormalizations 
are basic requirements following from the commutation 
relations of quantum spin operators, and are further dis- 
cussed later. 

As all calculations are done with RPA-lcvel magnon 
energies and spectral functions, this provides a RPA-levcl 
description of the magnetic phase diagram. Features such 
as continuous spin-rotational symmetry and transverse 
spin fluctuations are included, resulting in substantial 
improvement over the HF theory. Proper limiting results 
are therefore obtained as d ^ 2 (Tn — > 0, for all J7), and 
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d ^ oo (Tn ^ T^^, for U/t oo). 

We first discuss the approach for determining the Neel 
temperature within the spin picture used here in terms 
of the local spin operators Si. In the anisotropic broken- 
symmetry state for T < Tn, we have (SfS'f) > {Sf)/3. 
However, with increasing temperature the anisotropy 
progressively decreases, and eventually vanishes at the 
Neel temperature. Therefore we determine Tn from the 
following isotropy requirement, 

{SrS+) = {S+Sr) = {2/3){Sf) as T ^ . (1) 



I. STRONG COUPLING LIMIT 

We self consistently obtain the enhancement in trans- 
verse spin fluctuation {S^ S^) due to thermal excitation 
of magnons. The amplitude of the magnon propagator 
and the magnon energy are both multiplicatively renor- 
malized by (5'^), the sublattice magnetization. This fol- 
lows from a finite-temperature extension of the one-loop 
correction to the magnon propagator, obtained earlier 
in the strong correlation limit. |T3] Alternatively, in the 
equivalent spin-5' Heisenberg model, this can be seen 
within the equation of motion approach, which has 
been applied recently to the S = 1/2 case of layered an- 
tiferromagnets. ||l5| Considering first the simpler case of 
strong coupling, we therefore have for the transverse spin 
fluctuation, 
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where P — 1/T with fee set to 1, the magnon energy 



For S = 1/2 this yields Tn = 3.95t'^/U, in good agree- 
ment with the result 3.83t^ /U from high-temperature se- 
ries expansions. Q Furthermore, as the Goldstone mode 
is built-in, the spin fluctuation theory also yields the ap- 
propriate T^-falloff of the sublattice magnetization at low 
temperatures (T << Jn), and the vanishing of Tn in two 
dimensions, where fsp diverges. 

Interestingly, for hypercubic lattices this spin- 
fluctuation correction exactly matches with the Onsager- 
correction to the mean- field result. Here one obtains 
the reaction-field correction factor /orf = J2q ^/(l+Tg): 
which is identically equal to the spin-fluctuation correc- 
tion factor /sF for hypercubic lattices. It is thus quite 
a coincidence that while the (Onsager-corrected) MFT 
and SFT yield very different falloff of staggered magne- 
tization with temperature, identical Neel temperature is 
obtained in both theories for all S. 

The Onsager correction has recently been utilized to 
discuss the finite-temperature magnetism in the Hubbard 
model. Here the low-energy excitations of the Hub- 
bard model are approximately mapped to those of an 
effective 5 = 1/2 Heisenberg model with [/-dependent, 
extended-range spin couplings. The subsequent use of 
mean field theory, however, does not properly account 
for the continuous spin-rotational symmetry of the sys- 
tem, resulting in an expression for Tn which is insensitive 
to whether the effective spin model is of the Heisenberg 
or Ising type. As a vanishing Neel temperature would 
be obtained even for the Ising AF on a square lattice, 
the divergence of the reaction-field correction in two di- 
mensions should therefore be seen as a breakdown of the 
ORF theory in two dimensions, rather than a source of 
the vanishing of Tn. Consideration of the gapless spin- 
fiuctuation modes is essential for this. 



II. GENERAL CASE 



zJ (S"^} Jl - 72, and the two terms are respec- a. Evaluation of transverse spin fluctuations (T = 0) 



tively the thermal and quantum (zero-point) contribu- 
tions. (S/^S^) is given by a similar equation except 
with [(1 — 7g)^^/2 + 1] for the quantum contribution. 
The renormalization of the magnon propagator by the 
magnitude (25^) in Eq. (2) ensures that the commu- 
tation property of spin angular momentum operators 
= 2S^ is obeyed by the expectation value in 
the spin state. As (S^) vanishes in the limit T T-^ , 
from the isotropy requirement we obtain the Neel tem- 
perature. 



Tn = zJ 



S{S+1) 



SF 



(3) 



where the spin-fluctuation correction factor /sf = 
Y.q 1/(1 - Iq) in terms of 7, = Y.t=i cosg^/d. For the 
simple cubic lattice /sf = 1.517. As the mean-field result 
is T^'^ = zJS{S+l)/3, for the spin-S* Heisenberg model, 
spin fluctuations reduce Tn to nearly 70% of its MF value. 



In order to extend these results to arbitrary U/t, 
we need the magnon mode energies and spectral func- 
tions, and we now describe a new method for evaluat- 
ing transverse spin fluctuations at the RPA level, from 
which both the Neel temperature and the quantum spin- 
fluctuation correction to sublattice magnetization are 
obtained. In addition to the pure antiferromagnets in 
two and three dimensions, this approach is applicable 
to a variety of other systems such as AF with disor- 
der, defects, vacancies etc., and provides a simple way 
to obtain first-order spin-fluctuation corrections about 
the HF-level broken-symmetry state. The method is 
based on a convenient way to perform the frequency in- 
tegral to obtain spin correlations from spin propagators, 
and we illustrate it here for transverse spin correlations. 
We write the time-ordered transverse spin propagator 
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(*G|T[S',r(OS'+(t')]|*G) at the RPA level in frequency 
space as, 



|(/)„(a;))((^„(u;)| 
(4) 



where (^„(a;),A„(w) are the eigensolutions of the x'^(w) 
matrix. Here x°(a-') is the zeroth-order, antiparallel- 
spin particle-hole propagator, evaluated in a suitable 
basis in the self-consistent, broken-symmetry state at 
the HF level. In the site representation for instance, 
[X°{u;)],j ^ i ^{duj' /2tt)G\^{uj')G^.^{lj' -uj) in terms of the 
HF-level Green's functions. Spin correlations are then 
obtained from, 



{s;{t)si{t')) 



cLu 
2^ 



where the collective mode energies LUn are obtained from 
1 — UXn{LOn) — 0, and A„((jj) has been Taylor-expanded 
as A„(w) « Xnit^n) + (w — uJn)idXn / dLL>)uj,^ near the mode 
energies to obtain the residues. For convergence, the re- 
tarded (advanced) part of the time-ordered propagator 
X ^, having pole below (above) the real-w axis, is to be 
taken for t' < t it' > t). The frequency integral is con- 
veniently replaced by an appropriate contour integral in 
the lower or upper half-plane in the complcx-w space for 
these two cases, respectively, which results in Eq. (5). 

For the pure AF case, it is convenient to use the two- 
sublattice representation due to translational symmetry, 
and we work with the 2x2 matrix [^''(qa;)], which is 
given in terms of eigensolutions of the HF Hamiltonian, 
p8| and numerically evaluated using a momentum grid 
with Afc — 0.1. Equal-time, same-site transverse spin 
correlations are then obtained from Eq. (5) by summing 
over the different q modes. We consider t' , so that 
the retarded part is used, with positive mode energies. 
From spin-sublattice symmetry, correlations on A and B 
sublattice sites are related via (5+5*^)^ = {S^S^)b- 



{S') = S 



1 f {s+s-) + {s-s+) 

S\ 2S 



1/2 



(6) 



To order 1/25, this yields the correction to the sublattice 
magnetization of {2S)-^ E^Kl " 7,)"^^^ - 1] = 0-156 in 
three dimensions for S = 1/2. This same result at one- 
loop level was earlier obtained from a different analysis in 
terms of the electronic spectral-weight transfer, JO] and 
is in exact agreement with the SWT result. |20|j21|| 

For arbitrary U, as discussed above Eq. (1), S = 
{S^)nF plays the role of the effective spin quantum num- 
ber. The sublattice magnetization m is therefore ob- 
tained from m = mfjF ~ ^™SF: where the HF sublat- 
tice magnetization is obtained from the self-consistency 
condition, [ p^ and the spin-fluctuation correction Smsp 
at the one-loop level is obtained from Eq. (6) with 
S = {S')hf, 



6msF = 



{s+s- + S-S+) 



RPA 



{S+s- -S-S+)npA 



1. 



(7) 



The [/-dependence of the sublattice magnetization m{U) 
is shown in Fig. 1, showing that it properly interpo- 
lates between the weak and strong coupling limits, ap- 
proaching the SWT result, m{oo) = toswt = 0.844 as 
U/t oo. As expected, the spin-fluctuation correction 
(5tosf increases with interaction strength. We note that 
the two \ow-U datapoints are subject to the same finite- 
size numerical error which affects the transverse spin cor- 
relations, as mentioned earlier. 

Earlier studies of m{U) have employed self-energy cor- 
rection {d = 2, mix,) ~ 0.7 ^ mswx), (H the QMC 
method (d — 2), p3] functional-integral schemes (d = 3, 
m(oo) — 1.0), and recently an approximate mapping 
to an extended-range QHAF with effective U -dependent 
spin couplings obtained from the static transverse spin 
propagator at the RPA level {d = 3, m{oo) = mswx)- 
In the last work, use of the static propagator im- 
plies that dynamical effects are not fully included, which 
arc in principle important in weak and intermediate cou- 
pling regimes when the magnon energy is not negligible 
compared to the energy gap. 



B. Quantum correction to sublattice magnetization 



C. Finite temperature 



From the commutation relation [5+,S' ] ~ 25*^, the 
difference {S+S~ - S'-S'+)rpa should yield (2S'^)hf, 
which is indeed confirmed as shown in Fig. 1. The de- 
viation at small U is presumably because of the finite 
grid size used in the calculations. The sum {S+S~ + 
S~ S+)^pA yields a measure of transverse spin fiuc- 
tuations about the HF state, and for spin S in the 
strong coupling limit, one obtains {S+ S^ + S^ S+)-rpa = 

(2^) 1/^1 -7g'- Using the identity, (S^S^) = 

S{S+1) - (S+S- +S-S+)/2, the sublattice magnetiza- 
tion m = (25^) is then obtained from, [pO[ 



To extend these calculations to finite temperature, we 
write the diagonal elements Xu^i^) and X^~(w) in terms 
of the RPA- level spectral functions and mode energies as. 



-Ar, 



B„ 



(8) 



and a similar equation for xti~{^) with Aq and Bq in- 
terchanged, the symmetry following from the identity 
X ^(ijj) = (~^)- In the strong coupling limit, the 
spectral functions have the familiar form: 
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FIG. 1. The sublattice magnetization m vs. U (diamonds), 
along with the HF results from (i) the self-consistency condi- 
tion (dashed), and (ii) (2S")hf = (S+S" - S'"S'+)rpa (plus). 
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l]/2 and i3, = [(l-72) 



-1/2. 



l]/2. 



The transverse spin correlations at finite temperature are 
then obtained conventionally from these propagators by 
using the Poisson summation formulas which allow re- 
placing the sum over Matsubara frequencies by a contour 
integral in complex frequency space, and by picking con- 
tributions from both poles in the propagators, traversed 
in the clockwise direction. |23| We obtain. 



-A„ 



B„ 



(9) 



Both Aq and Bq are determined from the retarded 
part of X ^(oj) having positive mode energies by us- 
ing the spin-sublattice symmetry x~aaX~^^ = xJia(^) ~ 
Xj3^(w). Thus the spectral functions are obtained from 
the magnon amplitudes on A and B sublattices using. 



' U^{dXq/dLu) 



B\2 



Bq = 



(10) 



In analogy with the strong-coupling result discussed 
above Eq. (2), and in keeping with the renormaliza- 
tion constraints of this RPA theory as discussed above 
Eq. (1), we now multiplicatively renormalize both the 
magnon energy ujq, and the magnon amplitudes Aq,Bq, 
by multiplying with the momentum-independent ratio 
(2S'^)/(2S'^)hf. This ensures that the commutation 
property [5"+, 5^] = 2Sz is satisfied by expectation val- 
ues in the spin state, as seen below. By separating 
the zero- and finite-temperature contributions in {S~ S~^) 
and {S^S^) from Eq. (9) and its counterpart, it is seen 
that ([5'+,S'-]} = J2qiBq - Aq) = (2S'^)hf at all tem- 
peratures. This HF result from use of the RPA-level 
transverse spin correlations is corrected by the above 
renormalization, and the exact result (25*2) is obtained 
instead. Furthermore, as in the strong coupling case, 
the renormalization of the magnon energy ensures that 



the spin stiffness vanishes at the Neel temperature where 
long-range order vanishes. However, this does not consti- 
tute an independent renormalization as, in fact, the two 
renormalizations of magnon energy and amplitudes are 
connected. The only renormalization in x^i^) (of the lu- 
term) which affects the magnon amplitudes also modifies 
the magnon energy in precisely the same manner. p3[ 



As T 
Eq.(9), 



and (S^) — > 0, we therefore obtain from 



{S-S+) = {S+S-)=TnY, 



Ar,+B„ 



(11) 



where Aq = ^,/(25")hf, Bq = B,/(25^)hf, and Cjq = 
ijjq/ {2Sz)bf are the scaled spectral functions and magnon 
energies. As for the strong coupling case, Tn is obtained 
from the isotropy requirement given in Eq. (1), with 
{S"^) — S{S -f 1) in terms of the effective spin quantum 
number S = {Sz)-hf = "^hf/2. With this the Neel tem- 
perature is then obtained from 



In 



E 



B„ 



(12) 



This result has the proper limiting behavior in both 
limits. In the weak coupling limit, the magnon energy 
ujq « 2A, the maximum magnon energy, for almost all 



q modes. Also Aq + Bq « 1, so that with mnF 



we obtain Tm 



0, 



the proper HF scale. 



And in the 

2M-1 



strong coupling limit, {Aq + Bq)/ujq = [3J(1 — 7^)] 
so that Eq. (3) is obtained, which is in good agreement 
with the result from high-temperature series expansion, 
as discussed earlier. 

Fig. 2 shows the variation of Tn with J7, as obtained 
from above equation, along with the HF results which is 
asymptotically approached from below in weak coupling. 
In strong coupling the result closely approaches the high- 
temperature series-expansion result 3.83t^/C/, which is 
well below the MFT result 6t^/J7. For weak couphng 
HF-level finite-temperature effects need to be considered 
in Eq. (12) when the temperature is not negligible com- 
pared to the energy gap, and so it is the effective spin 
quantum number S = mHF/2 at temperature Tn which 
is used here. However, the magnon mode amplitudes and 
energies are obtained at T = 0. Determination and use 
of magnon spectral properties at temperature T = Tn is, 
of course, possible by iteration, but not very practical. 
We find that these HF-level finite-temperature correc- 
tions become negligible for U /t > 5 in the crucial strong 
correlation limit. On the weak-coupling side of this, the 
energy gap and therefore the magnon energies ujq would 
be reduced somewhat at finite temperature, so that the 
results for Tn are a slight overestimate for U /t < 5. 
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FIG. 2. Neel temperature Tn vs. U /t (solid line), and us- 
ing the dominant- mode approximation (dashed line). Also 
shown are the strong-coupling asymptotes MF and HT, from 
mean-field theory (T^^ — Qt^ /U) and high-temperature se- 
ries expansion {T^'^ = 3.83t^ /U) respectively. 



D. Dominant-mode approximation 

We now consider a dominant-mode approximation 
which, in the strong coupUng limit, becomes equiva- 
lent to the dynamical mean-field theory. In the strong- 
coupling limit, where ujq = dJ^l — 7^ in d dimensions, 

the magnon density of states is strongly peaked at the up- 
per band edge at energy dJ, becoming a delta function in 
the limit d ^ 00. As the highest-energy modes are domi- 
nant, in this approximation we only consider these modes 
in the sum in Eq. (12). As charge fluctuations which 
have been ignored are absent in the strong coupling limit, 
and corrections to RPA decrease as 1/d, this dominant- 
mode approximation becomes exact in the limit of infi- 
nite dimensions. This approximation also becomes good 
in the weak-coupling limit, where the magnon spectrum 
is strongly peaked at the gap energy 2 A, and all fluctua- 
tions are weak. Since for the highest-energy modes with 



ujm, we have Aq 



Bq « 1, Eq. (12) simplifles to 



(13) 



which is also shown in Fig. 2. As LOm = 3 J in the strong- 
coupling limit U /t 00, this yields Tn = 3J/2 = 6t^/U, 
the mean- field result, which is also asymptotically ap- 
proached in the QMC ||] and DMFT H calculations. 
On the other hand, in the limit U/t — > 0, as LOm ap- 
proaches 2A, the energy gap, we obtain Tn ~ A, the 
proper weak-coupling scale. This analysis indicates that 
the essential features of the Tn vs. U behavior are con- 
tained in LOm, the energy scale of low-energy spin excita- 
tions associated with spontaneous symmetry-breaking of 
the continuous spin-rotational symmetry. A relationship 
between Tn and the low-energy scale was also obtained 
for the infinite-dimensional Bethe lattice. [E6| 



At finite temperature, the gap parameter A is de- 
termined from the self-consistency condition, 1/U = 
EJl/2£;fe)tanh(/3Sfc/2), where = T^M^ is the 
AF band energy. In the weak coupling limit, when the 
temperature is comparable to the interaction strength, 
the gap parameter, and therefore the maximum magnon 
energy, are sensitive to temperature. We therefore use 
tJm(C/, Tn) in Eq. (13), which is obtained from the 
magnon spectrum for a gap parameter A(Tn) at tem- 
perature Tn. The self-consistent determination of Tn(C/) 
is easily implemented in two steps — (i) obtaining Li;„i(A) 
and hence Tn(A) as a function of the gap parameter A, 
(ii) then using the self-consistency condition to obtain 
C/(A,Tn). As mentioned earlier, these HF-level finite- 
temperature corrections become negligible for \J ji > 5. 
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FIG. 3. Magnon energies coq vs. U at momenta {n,n /2,0) 
(solid) and (7r/2, 7r/2, 7r/2) (dashed), together with the Hub- 
bard gap 2 A (dotted). 

The [/-dependence of ojm at T = is shown in Fig. 
3, showing that it properly interpolates between the two 
limits. Similar qualitative nature of the U-dependence of 
the maximum magnon energy was found in two dimen- 
sions, ||2^ and infinite dimensions. [^,^ The complete 
spectrum of spin excitations was studied earlier for the 



d — 2 case. |27| The band-edge degeneracy at energy 
uj„i — 2 J in the Hcisenberg limit is removed for finite U, 
and the maximum was found to occur for q — (ivr, 0) and 
(0, ±7r). For d = 3 we find that for most of the U/t range 
the maximum occurs at momenta q = (±7r, ±7r/2, 0) 
etc. where all the cosg^ terms are completely different. 
For small U (less than about half bandwidth B/2) the 
magnon energy at (7r,7r/2, 0) is actually not the maxi- 
mum, as shown in Fig. 3. However, the difference is 
marginal, and we have neglected this detail. As the de- 
generacy of g-points of the type (tt, 7r/2, 0) is much higher, 
yielding a correspondingly higher DOS, it is reasonable 
to use this magnon energy as the relevant energy scale 
uj,n- This also implies that the structure of the magnon 
density of states undergoes a crossover (at C//t « 6) , and 
for small U the peak in magnon density of states changes 
position from the upper band edge to within the band. 
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FIG. 4. The mode-averaged ratio i? = ^ (um/i-^q) vs. U , 
showing the cusp at U/t ~ 6, which arises from the red-shift 
at lower-f/ of the magnon DOS-peak energy relative to the 
maximum magnon energy ujm- 



properly between the weak and strong coupling limits. 
This analysis and the new method to evaluate transverse 
spin fluctuations should also prove valuable to other re- 
lated systems such as antiferromagnets with disorder, de- 
fects and vacancies, where spin-fluctuation studies can be 
extended to strong disorder and strong defect/ vacancy 
concentrations. 
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This crossover has been seen in a recent study of the 
[/-dependence of the magnon density of states in d = 3. 




A convenient way to observe this crossover is through 
the [/-dependence of the mode-averaged ratio R = 
J^qi^m/tOq). Hcrc LOm is the maximum magnon energy 
(band edge), obtained by sampling through the full spec- 
trum. R approaches 1 in the limit when all magnon 
modes have energy approaching ujm- The [/-dependence 
of R is shown in Fig. 4. With decreasing U from the 
strong-coupling limit, there is at first a marginal increase 
in R, which then rapidly starts decreasing a,t U ~ B. 
This decrease in R indicates a reduction in the magnon 
dispersion, so that relatively more modes have energy 
closer to the maximum magnon energy u>m- This trend 
is, however, sharply cutoff at [/ « 5/2 = 6t, as clearly 
seen from the cusp in Fig. 4. The sharp enhancement in 
R with decreasing U/t is due to the onset of a red-shift 
of the DOS-peak energy (where the magnon DOS peaks) 
relative to the band edge, uJm- Because of this red-shift 
most magnon modes have energy less than ujm, yielding 
an enhanced contribution to R. The rapid decrease in R 
towards unity with further decrease in U indicates that 
almost all magnon modes have energy nearly equal to 
UJm which, as mentioned earlier, approaches the charge 
gap 2A in weak coupling. The role of the approaching 
Stoner band in compressing the magnon spectrum has 
been noted earlier. 

In conclusion, this analysis in terms of transverse spin 
fluctuations at the RPA level provides a good account of 
the magnetic phase boundary for the three-dimensional, 
half-filled Hubbard model, particularly in the crucial 
strong coupling limit. The quantum correction to sub- 
lattice magnetization as well as the Neel temperature are 
obtained in the whole U/t range, and both interpolate 
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